Gauge-invariant electromagnetic response of a chiral px + ipy superconductor 
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We present a gauge-invariant theory of the electromagnetic response of a chiral p^ + ipy super- 
conductor in the clean limit. Due to the spontaneously broken time-reversal symmetry, the effective 
action of the system contains an anomalous term not present in conventional superconductors. As a 
result, the electromagnetic charge and current responses contain anomalous terms, which explicitly 
depend on the chirality of the superconducting order parameter. These terms lead to a number of 
unusual effects, such as coupling of the transverse currents to the collective plasma oscillations and 
a possibility of inducing the charge density by the magnetic field perpendicular to the conducting 
planes. We calculate the antisymmetric part of the conductivity tensor (the intrinsic Hall conduc- 
tivity) and show that it depends on the wave vector of the electromagnetic field. We also show that 
the Mermin-Muzikar magnetization current and the Hall conductivity are strongly suppressed at 
high frequencies. Finally, we discuss the implications of the theory to the experiments in Sr2Ru04. 

PACS numbers: 74.70.Pq, 78.20.Ls, 74.25.Nf, 73.43.Cd 



I. INTRODUCTION 

Unconventional superconductors with spontaneously 
broken time-reversal-symmetry (TRS) recently attracted 
significant interest [1, 2]. The idea of a superconduct- 
ing pairing violating additional symmetries of the normal 
phase (on top of the gauge symmetry) is intriguing, and 
there has been a lot of effort to find materials exhibit- 
ing such pairings. Considerable evidence indicates that 
Sr2Ru04 is an unconventional superconductor with bro- 
ken TRS [3-6]. The most convincing indication is the 
recent observation of the polar Kerr effect (PKE) in the 
superconducting state of Sr2Ru04 by Xia et. al. [3]. In 
this experiment, a linearly polarized light, incident along 
the z direction perpendicular to the conducting planes 
of Sr2Ru04, experiences normal reflection [7]. R was 
found in Ref. [3] that the polarization plane of the re- 
flected light is rotated relative to the polarization plane 
of the incident light by the Kerr angle Ok 65 nrad. 
The effect appears below the superconducting transition 
temperature Tc — 1.5 K. The Kerr rotation, which may 
be clockwise or counterclockwise, develops in the absence 
of an external magnetic field and is a clear signature of 
the spontaneous TRS breaking in the superconducting 
state. The experiment [3] used the Sagnac interferom- 
eter, where two counterpropagating laser beams retrace 
their paths, so that all effects, other than the TRS break- 
ing in the sample, cancel out exactly. Although previous 
muon-spin- relaxation measurements [5, 6] gave an indi- 
rect evidence that the TRS is broken in Sr2Ru04, the 
PKE experiment [3] provides much stronger evidence for 
this remarkable effect. Additional indication for the TRS 
breaking in Sr2Ru04 comes from the Josephson junction 
experiment in the presence of a magnetic field [4] , which 



was interpreted as the evidence for existence of domains 
with opposite chiralities. On the other hand, the scan- 
ning SQUID and Hall probe experiments [8, 9], designed 
to search for domains with opposite chiralities at the sur- 
face of Sr2Ru04, did not find any evidence for the TRS 
breaking. These results show that macroscopic manifes- 
tations of the microscopic TRS breaking are not fully 
understood and require further theoretical investigation. 
In this paper, we study the electromagnetic properties of 
a chiral superconductor with the -\- ipy pairing in the 
clean limit. 

Sr2Ru04 is a layered perovskite material consisting of 
weakly coupled two-dimensional (2D) metallic sheets par- 
allel to the {x,y) plane [10, 11]. It was proposed theoret- 
ically that the superconducting pairing in this material 
is spin-triplet [12] and has the chiral Px + ipy or p^ ~ ipy 
symmetry [13]. In this state, Cooper pairs have the or- 
bital angular momentum = -\-h or = —h normal to 
the layers. Such an order parameter breaks the TRS and 
is analogous to the 2D superfiuid ^He-A [14]. It should 
be emphasized that the questions of the spin symmetry 
(singlet vs triplet) and the orbital symmetry (chiral vs 
nonchiral) are separate issues. It is possible to construct 
chiral order parameters for both triplet and singlet pair- 
ing [15]. There is substantial experimental evidence in 
favor of the spin-triplet and odd orbital symmetry of pair- 
ing in Sr2Ru04 [10], which includes measurements of the 
spin susceptibility [16-18] and the Josephson effect [19]. 
However, there are also alternative interpretations [15] in 
terms of singlet pairing. We study the electromagnetic 
response of quasi-2D (Q2D) chiral superconductors, and 
our results are applicable (with minor modifications) for 
either spin symmetry. We do not pay special attention 
to the existence of nodal lines in the order parameter of 
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Sr2Ru04 (sec, e.g., Ref. [20] for an interpretation of tun- 
neling measurements). The nodal lines do not affect the 
chiral response qualitatively, so we concentrate on the 
simplest case of the Px + ipy pairing. 

Although the Px + ipy superconducting pairing breaks 
the TRS and, in principle, permits a non-zero Kerr angle 
6k, an explicit theoretical calculation of 9k is challeng- 
ing. A textbook formula [21] expresses 9k in terms of the 
ac Hall conductivity axyi^) at the optical frequency lj. A 
calculation of the intrinsic Oxy{uj) for a chiral supercon- 
ductor in the absence of an external magnetic field turns 
out to be quite nontrivial. It is customary for theoretical 
calculations to use the gauge where the scalar potential 
Aq is set to zero and only the vector potential A is con- 
sidered. In this gauge, the calculations [22, 23] show that 
there are no chiral terms in the single-particle response of 
a + iPy superconductor. Using this gauge and taking 
into account particle-hole asymmetry, Ref. [24] found a 
small chiral response from the collective flapping modes. 

However, when calculations are performed in a gen- 
eral gauge, they do produce a non-trivial Chern-Simons- 
type (CS) term in the effective action of a Q2D px + ipy 
superconductor, which breaks the TRS: 

5cs = ey dtd^r {Ao + dt^/2e) (dyAx - dxAy)/c. (1) 

Here, $ is the phase of the superconducting order pa- 
rameter; e and c are the electron charge and the speed 
of light, respectively. Equation (1) was first derived in 
Ref. [14] at T = 0, and the coefficient Q was found to be 
6 = ±e'^/2hd, where h and d are the Planck constant and 
the distance between the layers. The signs ± correspond 
to the Px ± ipy pairing. The term (1) was then studied 
in Refs. [25-30], either at T = or at T close to % using 
the Ginzburg-Landau expansion. In a recent paper [31], 
Sqs was obtained for a finite frequency lu and arbitrary 
temperature T, in which case Eq. (1) should be written 
in the Fourier representation with the coefficient <d{u>) 
under the integral. Ref. [31] found that 6(w) decreases 
as (Ao/?iw)^ at high frequencies huj 3> Aq, where Aq is 
the superconducting gap [see Eqs. (71)-(74)]. 

Equation (1) is the only term in the effective action of a 
Px + ipy superconductor that breaks the TRS. Under the 
time-reversal operation, the variables in Eq. (1) trans- 
form as A — > —A, Aq — > dt — > —dt, and $ — > — 
so Scs changes sign. The coefficient 9 in Eq. (1) ex- 
plicitly depends on the chirality of the order parameter 
Px ± ipy and changes sign when the time-reversal opera- 
tion is applied to the superconducting pairing itself. The 
action ^cs is similar to the standard Chern-Simons term, 
which has the structure ri^jj\A^d,jAx, with ri^,yx being 
the (2-|-l)D antisymmetric tensor and the indices /U, u, 
and A taking the values t, x, and y. Compared with the 
full, gauge- invariant Chern-Simons action, Eq. (1) misses 
the term AxdtAy — AydtAx- In Eq. (1), the gauge invari- 
ance is ensured by the presence of the superconducting 
phase <I>, which changes upon a gauge transformation to 
compensate for the change of Aq [14] . 



Taking a variation of the effective action with respect 
to A, one obtains the electric current j = —cdS/SA. For 
a full Chern-Simons term, this would give the usual Hall 
effect jx = UxyEy. However, a variation of Eq. (1) gives 
the following (anomalous) current 

= -e i X V(Ao + dt^/2e), (2) 

where f is a unit vector perpendicular to the conducting 
planes. It was shown in Ref. [29] that the current (2) can 
be expressed as a Mermin-Muzikar current [32] 



where p is the electron charge density, and rUe is the 
electron mass. Equation (3) can be understood as the 
magnetization current J = V x M originating from the 
magnetization M = —{fip/4me)z associated with the 
angular momentum L = hz oi each Copper pair [29, 32]. 
While Eq. (3) is valid at low frequencies, our calculations 
show that at high frequencies it is suppressed by a factor 
{Ao/hu))"^, see also Ref. [33]. Moreover, we show that, in 
addition to the magnetization current (3), there is also 
an anomalous electric polarization current, that gives a 
comparable contribution even at low frequencies. 

Equation (2) is similar to the standard expression for 
the Hall conductivity, but its right-hand side does not 
contain the complete electric field E = —'VAq — dtA/c 
and includes the superconducting phase Equation (2) 
can be expressed in terms of E as 

a[V$- (2e/c)A]\ 

The second term in the brackets of Eq. (4) is propor- 
tional to the acceleration of the London supercurrent 
js = (ps/2me)[V<I> — {2e/c)A], where ps is the super- 
fluid charge density. It was argued in Ref. [31] that this 
term in Eq. (4) becomes ineffective at high frequencies 
and may be neglected. Then, Eq. (4) reduces to the 
conventional Hall relation, and the coefficient G can be 
identified with the Hall conductivity a^y = 6 and used 
for a calculation of the Kerr angle 6k [31, 33]. 

However, within the two-fluid model of superconduc- 
tivity, one can argue that Newton's equation of motion 
for the supercurrent is (m/e) dijs = PsE, i.e., the su- 
percurrent is accelerated by the electric force. Then, the 
right-hand side of Eq. (4) vanishes, and the chiral Hall 
current j^"'^ is zero [34]. The reason for this cancellation 
is that the superconducting phase $ has its own dynam- 
ics and compensates the electromagnetic field in Eq. (2) 
[35-37] . In general, Eq. (2) should be supplemented with 
an equation of motion for and then $ should be elimi- 
nated, so that the current response is expressed in terms 
of the electromagnetic field only. In other words, one 
should derive the effective action S{A^, $) for apx + ipy 
superconductor as a function of the electromagnetic field 
A^ and the superconducting phase and then integrate 



out $ and obtain a new action S{Afj_) in terms of the 
electromagnetic field only. This procedure is well estab- 
lished for nonchiral superconductors [38-41], and it was 
implemented for chiral superfluids in Refs. [25, 26, 28]. 
However, the calculations in Refs. [25, 26, 28] were per- 
formed only in the limit of low frequencies. As a result, 
some terms in the effective action were neglected, and 
the frequency dependence of the coefficients in the ac- 
tion was not considered. In our paper, we perform a 
detailed derivation of the effective action by taking into 
account dynamics of $ and the internal Coulomb poten- 
tial. Our results are applicable for all frequencies and 
exhibit non-trivial frequency dependence. 

Our calculations show that the effective Hall conduc- 
tivity depends not only on frequency, but also on the 
wave vector and is proportional to the square of the wave 
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parallel to the layers: axy oc qf^. This 



conclusion agrees with Refs. [25, 26, 28]. By taking the 
limit of q ^ at w 7^ [42], we find that the Hall 
conductivity for a spatially homogeneous system van- 
ishes, i.e., the cancellation in Eq. (4) indeed takes place. 
This result is consistent with the general conclusion of 
Ref. [43], which argued that an electric field cannot pro- 
duce a sideways motion of the electron gas without an 
external magnetic field, no matter what the internal in- 
teraction between electrons producing the p-wave pairing 
is. This argument is based on the Galilean invariance and 
is applicable to an infinite, spatially homogeneous, clean 
system without boundaries and impurities. 

This conclusion does not contradict the results of the 
PKE experiment, because the setup used in Ref. [3] gen- 
erates spatial inhomogeneities within the (x, y) plane of 
Sr2Ru04. Indeed, as sketched in Fig. 1, Xia et. al. [3] 
used a tightly focused Gaussian laser beam of the trans- 
verse size / ~ 25 ^m, not an infinite uniform electromag- 
netic plane wave. (The beam size I is smaller than the 
sample size and the size of a domain with a given chi- 
rality.) In the Fourier representation, this means that 
the electromagnetic wave has non-zero in-plane Fourier 
components g|| of the order of l/l. Because 6k is pro- 
portional to axy oc qp we conclude that the Kerr angle 
should be inversely proportional to P, 



Ok oc l/l^ 



(5) 



The theoretical prediction (5) can be checked experimen- 
tally. It is commonly assumed in literature that 9k de- 
pends only on lo and on the properties of a material [21]. 
However, Eq. (5) shows the Kerr angle for a chiral su- 
perconductor also depends on the geometrical size of the 
laser spot. 

Equation (3) shows that the chiral current oi & px + 
ipy superconductor directly couples to charge response 
and plasma collective modes. The Chern-Simons-type 
term (1) couples the longitudinal and transverse elec- 
tromagnetic fields and leads to a number of interesting 
effects, that are not present in nonchiral superconduc- 
tors. By taking a variation of Eq. (1) with respect to 
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FIG. 1: (color online). Schematic picture of the experimental 
setup in Ref. [3]. The incoming laser beam is incident nor- 
mally to the conducting planes of the layered superconductor 
Sr2Ru04. The beam is focused by the lens to a spot of the 
diameter I at the surface of the sample. The solid (red) line 
shows the Gaussian profile of the polarized electric field Ex 
in the beam (Ref. [7]). The interlayer distance in Sr2Ru04 is 
d. 



Aq, we find an anomalous charge response to the trans- 
verse electromagnetic field: ^p^"^ (x B^, i.e., the charge 
is induced by the magnetic field Bz along the z direc- 
tion [27, 29]. In Sec. VII B, we propose an experimen- 
tal setup to directly verify this relationship. We also 
calculate how the relationship between Jp'^"' and B^ is 
modified at high frequencies. By the continuity relation, 
the induced charge (Jp'"' produces an electric polariza- 
tion current, which gives an additional contribution to 
the Hall effect of the same order as the magnetization 
current. 



The paper is organized as follows. In Sec. II, we rig- 
orously derive a general expression for the effective ac- 
tion of a Q2D px + ipy superconductor and obtain the 
gauge-invariant electromagnetic response. In Sees. III- 
V, we discuss the collective modes and the conventional 
and anomalous (chiral) electromagnetic responses. In 
Sec. VI, we obtain the symmetric and antisymmetric 
parts of the conductivity tensor of a chiral superconduc- 
tor. In Sec. VII, we discuss the relationship of our results 
with the experimental studies of the chiral response of 
Sr2Ru04. Finally, we summarize the results in Sec. VIII. 
Some technical details are relegated to the Appendixes. 
In particular, a simplified alternative derivation of the 
effective action is given in Appendix C. 
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II. EFFECTIVE ACTION FOR A CHIRAL 
QUASI-TWO-DIMENSIONAL + ipy 
SUPERCONDUCTOR 

A. Triplet + ipy pairing 

First, we briefly summarize basic information about 
the triplet p-wave pairing [44]. The Cooper pairing 
between electrons is described by the pairing poten- 
tial Aa0{r,r') oc {ipa{r)'tlj0{r')). Here, tpair) is the 
electron annihilation operator at the point r with the 
spin projection a =t,i- For a uniform, translationally- 
invariant system, the pairing potential depends only on 
the relative distance r — r', so one can perform the 
Fourier transform and use the momentum representa- 
tion Aaj3{p) oc {'ipa{p)^0{—p))- The pairing tensor 
^affip) = A{p) d{p) ■ ct2 ??7/3 can be written in terms of 
the antisymmetric metric tensor rj^p and the Pauli ma- 
trices cr^, where the unit vector d{p) characterizes the 
spin polarization of the triplet state. The prefactor A(p) 
is a momentum-dependent pairing amplitude. 

For Sr2Ru04, we consider the case where the vec- 
tor d{p) = z has the uniform, momentum-independent 
orientation, which represents pairing between electrons 
with the opposite spins {ipi{p)ipi{—p)) . (It can be trans- 
formed into pairing with parallel spins by changing the 
spin quantization axis.) For the orbital symmetry, we 
consider the chiral pairing potential A(p) = Aq{px ± 
ipy)/pF, where pp is the Fermi momentum, and Aq is 
the superconducting gap. This order parameter corre- 
sponds to a vortex in the momentum space, because the 
phase of A(p) changes by ±27r when p goes around the 
Fermi surface. It is instructive to write the pairing po- 
tential in the form 



A{P) = A^P^ + iAy 



(6) 



and set A^ = ±Ay = Aq/pf only at the end of the 
calculations. The sign of the product 

Sxy = sign(Aa,Ay) (7) 

reflects the sign of the order-parameter chirality. 

B. Theoretical model 

Our goal is to derive an effective action for a chiral 
superconductor in a weak electromagnetic field. This ap- 
proach is equivalent to a linear response calculation [45- 
49]. We use the Greek indices, e.g., /x and v, to denote 
the space-time components of tensors in the Minkowski 
notation and the Roman indices, e.g., k and I, for the 
space components. To simplify intermediate steps of 
calculations, we set the Planck constant h, the Boltz- 
mann constant ks, and the speed of light c to unity: 
h = ks = c = 1. The constants % and ks can be easily 
restored in the final equations by dimensionality. The 
speed of light can be restored by noting that it appears 



only in the combination A/c with the vector potential A. 
To simplify presentation, we first study a purely 2D case, 
which corresponds to just one metallic layer in the (x, y) 
plane, and then generalize the calculation to the case of 
many coupled parallel metallic layers, as appropriate for 
Sr2Ru04. 

The Hamiltonian of interacting electrons, subject to an 
external electromagnetic field A^^, is given by 



H 



2m« 



+ j dVeAo(r)(5n(r) (8) 

- j d'^r d'^r' g{r - r')i}\{r)'fl}\{r')i}i{r')'il)-^{r) 
g2 r 

+ — (frcfr'6n{r)V{r-r')6n{r'). 

Here, rrif, and fj, are the electron mass and the chemical 
potential, g{r — r') is an anisotropic interaction potential 
leading to a p-wave pairing, and V{r — r') — l/|r — 
r'l is the Coulomb interaction potential. The density 
fluctuation operator 5n{r) reads 



6n{r) = tl;l{r)il)„{r) - uq, 



(9) 



with no being the 2D background charge density. The 
summation over repeated indices is assumed everywhere. 

The starting point for our calculation is the partition 
function Z, which can be expressed as a path integral 
over the anticommuting fermionic fields ip and ifj^ . We 
use the Hubbard-Stratonovich transformation to decou- 
ple the fermion interaction terms in Eq. (8) by introduc- 
ing additional integrals over auxiliary bosonic fields A 
and f [38-40, 50, 51]. The complex field A(r,r',r) is 
the superconducting pairing potential, and ip{r, r) is the 
internal electric potential produced by electrons. As a 
result, the partition function, written in the imaginary 
time T, reads [52] 



Z = j VA*VAVip e"^"- j Vtp^Vtp i 



(10) 



where the bosonic action is 



S... = ieJdrd^rHr,r)^Aoir,r)]no 



(11) 



g{r-r') 
and the electronic action is 

' [p - eA{r, r)]2 



2V{r-r') 



■I 



dTd'^rtp'^{r, t) 



2me 



+/<i..^rV.<(...K3.-i.|.(r..)+^»(r..)|)*,(r,.) 
-jdrd^rdPr' A(r,r',r)^|(r,r)^|(r',T)-FH.c. .(12) 



5 



The superconducting pairing potential in Eq. (12) can be 
written as a function of the relative coordinate r — r' and 
the center-of-mass coordinate R = (r + r')/2. For a uni- 
form system, in the absence of electromagnetic field, the 
equilibrium saddle-point configuration of Aq in the total 
action (10) does not depend on R and r, and is given 
by Eq. (6) in the Fourier representation with respect to 
r — r'. In the presence of an applied electromagnetic 
field, the saddle-point value of the complex supercon- 
ducting pairing potential is shifted and can be written in 
the adiabatic Born-Oppenheimer approximation [40, 53] 
as 

A{p, R, t) = Pl±iP^ [Ao + 5A{R, t)\ e^*(«'-) . (13) 
Pf 

We assume here that the applied electromagnetic field is 
weak, and varies slowly in space and time, i.e., ^/A <C 1 
and LO <C Wcutoflf < £f- Here, A and to are the wavelength 
and frequency of the electromagnetic field, ^ = vf/Aq 
is the superconducting coherence length, Wcuioff is the 
cutoff frequency of the interaction responsible for super- 
conducting pairing, and ef is the Fermi energy. The field 
shifts the amplitude of the pairing potential by 5 A and 
gives a space-time dependence to the phase $ of the or- 
der parameter. As in the s-wave superconductors, the 
amplitude variations SA are massive [47], and their con- 
tribution to the linear electromagnetic response is small 
in the parameter (Ao/e_F)^ <C 1 [46]. Therefore, we ne- 
glect the amplitude fluctuations in the rest of the paper 
and only consider the dynamics of the phase The p- 
wave order parameter also has other modes, such as the 
clapping modes, due to its internal orbital structure [24] . 
However, these modes are also massive, and we do not 
consider them. 

The phase $ of the order parameter is essential for 
ensuring the gauge invariance of the theory. By perform- 
ing a unitary transformation of the fermion operators 
^(j(r,r) e'*(^''^)/^Vcr(^i '^)i one can compensate the 
phase of the order parameter [40]. As a result, the elec- 
tromagnetic field in Eq. (12) is replaced by the gauge- 
invariant combinations of electromagnetic potentials and 
the superconducting phase: 



Here, ^ represents the integration over momenta as 
well as the summation over the Matsubara frequencies 
[54]. The terms G^^, Fi, and F2 in the electron ac- 
tion (15) contain the zeroth, first, and second powers of 
the electromagnetic fields and (p, 

G~^{p)=iuJra - £,ph -pxA^Ti +PyAyT2, (16) 

Fi(p, q) = -ie[Ao{q) + ^{q)]h - ev{p + q/2) ■ A{qll7) 



2 dpkdpi 



Ak{q)Ai{-q)T3. 



(18) 



Here, the wave vector q of the electromagnetic field and 
its bosonic Matsubara frequency f2„ are combined into 
q = (ir2„,q). The electron dispersion and velocity are 
= p^/2me — /X and v{p) = p/m^. Equations (16) (18) 
are written in terms of the Pauli matrices f acting on the 
electron spinor tj). 



C. Integrating out fermions 

Substituting the action (15) into Eq. (10) and taking 
the functional integral over ip, we obtain the electron 
contribution Sf^^ to the effective action up to the second 
order in the electromagnetic fields and (p. 



Sfs{A^, if) = -TVln(G-i + Fi + F2) 
- TrGFi - TrGFa + ^TrGFiGFi. 



-TrlnG"^ 
(19) 



Here, Tr denotes both the matrix trace and the sum over 

the internal frequencies and momenta. 

The term TrlnG^^ in Eq. (19) and the term propor- 
tional to |Ap in Eq. (11) define the saddle point config- 
uration (6) of the unperturbed system. The functional 
integral over the amplitude of A is taken by expansion in 
the vicinity of this saddle-point configuration. Combin- 
ing the remaining terms in Eqs. (11) and (19), we obtain 
the effective phase-only action of the system 



/ 



I?$I?^exp[-5'eff(>,(^)], 



(20) 



Ao = Ao-dr^/2e, A = A-V$/2e. (14) S,s{A^,ip) = -J2 



In this way, one ensures that the gauge invariance is ful- 
filled at every step of the calculation. 

After the phase transformation, the superconducting 
pairing potential in Eq. (12) has the equilibrium value 
given by Eq. (6). The electron Lagrangian (12) can 
be written in the momentum-frequency representation 
as a 2 X 2 Nambu matrix acting on the spinor xp{p) = 
[ip-^ {p, ujm) , i^li—p, oJm)], where ojm is the fermionic Mat- 
subara frequency, and p = {iu}m,Px,Py), 

Sei = J2^Hp)G-'{p)^{p) (15) 
p 

+ ^'>P^{P + 9)ri(p, q)^{p) + ^Hp)^2{p, q)ip{p). 



V{q) 



A{q)-M-Q) 



+-Tr[GFiGFi]. 



(21) 



P>9 



P,9 



Here, V{q) is the Fourier transform of the Coulomb po- 
tential, which in 2D is equal to V{q) = 27r/|q|. 

Now we need to evaluate the trace in Eq. (21). This 
amounts to calculation of one-loop Feynman diagrams 
using the fermionic Green's function G from Eq. (16), 

^ iWm + ^ph + PxAxfi - PyAyT2 , , 

and the vertex Fi (17), which represents interaction with 
the electromagnetic field. The final expression for the 
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effective action (21) can be written as 



(23) 



+iQok [io (<?) + P>{qMk i-q)+iQkoAk (q) [io (-g) + ^(-g)] • 

Here, Qoo, Ofei, and Qokiq) = Qko{-q) are the corre- 
sponding correlation (polarization) functions. 

The density-density polarization function Qoo = 

e2Tr[T3Gr3G] is 



00 = 



(fp 2e^ ■r-^iWm(za;m-|-zri„)-|-^_^+-A+A_ 



(27r)2 /3 



■E 



[(a;„ + n„)2+£;2][u;2^+i?2 



(24) 

Here, /3 = 1/T, e± = ep±q/2, A+A_ = A+A'+A+A", 

s± = ^e± + {Atr + {Atr with a±^, = (p.(,) ± 

9cc(j/)/2)A^(j^). 

The current- current correlation function Qki consists 
of the diamagnetic Q^;^ and paramagnetic Q^^i = 
e'^Ti[vkGviG] parts 



^fc; = ^ki, 



(25) 







(2). 
kl ' 



' (fp 2eVwi>r-i<^',u('-,» + '^^„j + i-4+ + A+A_ 



(27r)2 /? 



^E 



[(a;„ + Jl„)2+ii;2][a;^+£;2 



(26) 



with V = d^p/dp = p/rue being the electron velocity. 

The expressions (24), (25) and (26) for the density- 
density and current-current correlation functions are the 
same for chiral and nonchiral superconductors and con- 
tain nothing special. The only important difference be- 
tween chiral and nonchiral superconductors appears in 
the structure of the current-density correlation function 
Qko = e^Tr[wfcGT3G]. This difference plays the crucial 
role for the results of our paper. In a -|- ipy supercon- 
ductor, Qfco consists of the conventional, symmetric part 
Q^.q' and the anomalous, antisymmetric part Q^g'- 



QkO — Qk'o 



^kO 



^kO ' 



J (27r: 



cPp 2e'^Vk 
[ ll 

J (27r)2 /3 



(fp 2e^Wfe 



E 



E 



(27) 

[{uJm + ^nr + El][Ljl + Eiy 
(28) 

[iQxPy-iqyPx]^x^V 



[{Wm + ^nr + El][Ljl + Eiy 

(29) 

The two terms have the following symmetries 

where the operation q ^ —q means changing the signs 
of both frequency and momentum. Equation (30) fol- 

lows from an observation that Q^q is proportional to a 



product of the frequency and momentum components of 
q, whereas Q^q^ is a linear function of the momentum 
components q^ and qy only. One can also check using 
Eqs. (24), (25) and (26) that Qoo{-q) = Qoo{q) and 
Qkii-q) = Qki{q) are symmetric. 

Equation (29) shows that the anomalous charge- 
current correlation function Qj.q' explicitly depends on 
the chirality (7) of the order parameter, whereas the 
other correlation functions do not depend on it. We 
demonstrate in the rest of the paper that all of the non- 
trivial, chiral response of a Pa, -I- ipy superconductor origi- 
nates from the anomalous term (29). The density-current 
correlator Qko (27) is rarely discussed in textbooks and 
literature. Non-chiral superconductors have only the 
symmetric term Q^^^ (28), which couples to the longitudi- 
nal degrees of freedom (sec Sec. IV A), such as plasmon 
collective modes, but does not affect the transverse re- 
sponse, such as the London-Meissner current. Therefore, 
for the calculation of the transverse response, it is suffi- 
cient to consider a gauge with A-^ ^ and = 0. How- 
ever, in chiral superconductors, the antisymmetric term 
Q^^Q couples to the transverse response (sec Sec. VA) 
and controls the TRS-breaking response of the system. 
In the technical language, when calculating the trans- 
verse response of the chiral superconductors, one has to 
include vertex corrections in order to obtain correct re- 
sults. Without vertex corrections, calculations do not 
generate chiral terms in the electromagnetic response of 
&Px + iPy superconductor [22, 23]. 



D. Integrating out the internal Coulomb potential 

It is well-known that response functions of a charged 

superconductor, as opposed to a neutral superfluid such 
as ^He, are strongly modified by the Coulomb interaction. 
In our approach, integrating out the internal electric po- 
tential if is equivalent to taking into account the Coulomb 
interaction, because p mediates the electrostatic interac- 
tion between electrons. The effective action S'eff(A''), 
obtained after taking the functional integral over ip in 
Eq. (20), is defined as follows: 

Z = j V^V<pe-^'"^^"''^^ = j C^e-^'^f^"). (31) 

Substituting Eq. (23) into Eq. (31) and taking the Gaus- 
sian integral over p, we find 

5eflF(i'') = ]-^QQoiAo{q)iAo{-q) + QkiMq)M-q) 



+ 2QokiMq)Ak{-q), 



(32) 
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+ 




FIG. 2: Renormalization of the polarization functions (33) 
due to the Coulomb interaction V{q) in the RPA. The wavy 
line denotes V{q), which couples only to the charge vertices 
represented by the Pauli matrices T3. 



where the polarization functions are renormaUzed by the 
Coulomb interaction V{q) as follows: 



g, 



00 



l-V{q)Qoo{q)' 

QkoV{q) Qoi 



Q 



IkO 



kO 



l-V{q)Q 



00 



l-n9)Qoo(<z)' 

(33) 



The renormalized polarization functions (33) can be 
equivalently obtained through the resummation of the 
most diverging diagrams due to the Coulomb interaction 
in the random phase approximation (RPA) as shown in 
Fig. 2. 



F. Generalization to the quasi-two-dimensional 

case 

The derivation presented in the previous sections can 
be easily generalized to a Q2D system consisting of paral- 
lel superconducting layers separated by the distance d in 
the z direction. The layers are coupled via the interlayer 
electron tunneling amplitude tc- In the tight-binding ap- 
proximation, the electron dispersion and velocity are 



Px 



pI 



2me 



tccosip^d) - fi, v{p) 



dp ' 



(38) 



where (pxiPy) and arc the in-plane and out-of-planc 
momenta, respectively. We assume that the interlayer 
tunneling amplitude is much smaller than the in-plane 
Fermi energy: tc ^ ep- Then, the Fermi surface is 
a slightly warped cylinder extended in the Pz direc- 
tion. Equation (38) should be substituted into Eqs. (24)- 
(29), where the integrals should be taken over a three- 
dimensional (3D) momentum p, with p^ limited to the 
interval [— Tr/d, Tr/d]. 

The diamagnetic term Q^^-* in Eq. (25) should be re- 
placed by the expression following from Eq. (18) [40] 



_fp_ 

(27r)3 



1 _ ip tanh ( ^ 



E„ 



2T 



dpk dpi 



(39) 



E. Integrating out the superconducting phase 

We integrate out the superconducting phase $ to ob- 
tain the gauge- invariant effective action Sef[{A^^) that de- 
pends only on the electromagnetic field as follows: 



= / 



(34) 



where the field is defined by Eq. (14) and contains 
Substituting Eq. (32) into Eq. (34) and taking the 
Gaussian integral over we finally obtain the effective 
action [55] 



SMAn = ^^/^^.(gM''(9)A-(-9), 



Notice that the kernel (36) satisfies the identity 
q^K^, = K^.q'' = 0, 



(35) 



(36) 



(37) 



which ensures that the effective action (35) is gauge- 
invariant. Indeed, in the momentum representation, the 
gauge transformation is A^ — iq^4), where is an 

arbitrary scalar function. Substituting this expression 
into Eq. (35) and taking into account the identities (37), 
we see that the additional terms q^^cf) drop out and, thus, 
the effective action is gauge-invariant. 



where Ep ~ y Q, + The integrals (39) can easily 

be evaluated and give the following expression for the 
diamagnetic tensor of a Q2D superconductor: 



e^no ^ 



ikl 



nied 



nkU 




(40) 



Here, no = Pp/^n is the 2D electron density, and pp, 

vp, and ep = p|,/2me are the in-plane Fermi momen- 
tum, velocity, and energy. The dimensionless tensor n 
represents the anisotropy of the superfluid density. 

In a layered Q2D system, the renormalized polarization 
functions Q^j^ are given by Eq. (33) with the appropriate 
form of the Coulomb interaction [49] : 



V{q) 



2'Kd sinh(|q|||(i) 

|g||| cosh(|q|||ci) - cos{qzdy 



(41) 



where q|| = [q^, qy). Equation (41) reduces to the 3D ex- 
pression for the Coulomb potential Vm{q) at small mo- 
menta <C 1 and to the 2D Fourier transform of the 
Coulomb potential V2d(Q||) in the limit \q\\\d ^ 1 bs fol- 
lows: 



Vs-D{,q) = %, V2D(g||) = ?^. 

1 \q\\\ 



(42) 



In the long- wavelength limit (|q|d -C 1) considered in the 
rest of our paper, the appropriate form of the Coulomb 
potential is VzD{q). 
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G. Linear response 

By taking a variation of Eq. (35) with respect to A^, 
wc obtain the gauge-invariant electromagnetic response 

of the system, 



j^iq)=K^,Aq)A''iq). 



(43) 



Notice that the identities (37) ensure that the current 
(43) satisfies the continuity equation q^j^ = 0. 

To obtain physical results, we perform an analytical 
continuation from the Matsubara frequency to the real 
frequency u) + iS [52]. After the continuation, 

the vector q becomes g'^ = {u),q), and Eq. (43) gives 
a causal response to the external electromagnetic field 
= {Aq,A) at a finite temperature T. The kernel 
Knvisi) is defined via Eqs. (33) and (36) in terms of the 
one-loop response functions (24)- (29). Their analytical 
continuations to the real frequency after summation over 
the fermionic Matsubara frequencies are given in Ap- 
pendix A. 



III. COLLECTIVE MODES 



Here, the tensor n is defined in Eq. (40), and ujab and uic 
are the in-plane and out-of-plane plasma frequencies, 



2 A e^no 
w„h = 47r 



2 ^ 2.f!!l^^. (47) 



As shown in Appendix B, the collective mode (46) cor- 
responds to coupled plasma oscillations of the internal 
electric potential and the superconducting phase $. The 
plasma mode frequency (46) depends on the ratio of the 
in-plane q\\ and out-of-plane q^ momenta [57]. Plasma 
oscillations arc gapped in all directions, but arc strongly 
anisotropic due to the smallness of the interlayer tunnel- 
ing amplitude tc- The experimental values of the plasma 
frequencies in Sr2Ru04 are Wc = 0.32 eV and uJab = 4.5 
eV [58, 59]. 

If one formally sets ic = and uses the 2D expression 
for the Coulomb potential from Eq. (42), then Eq. (45) 
gives the plasmon mode dispersion for a single 2D layer 



' 2777106^ . 



me 



(48) 



Under certain conditions, an infinitesimal external 
electromagnetic field can induce a large current or den- 
sity response of the system, which indicates the existence 
of internal collective excitations in the superconductor. 
These resonances occur when the denominator in Eq. (36) 
goes to zero. Using Eq. (33), the denominator can be 
written as Q^i^q^q'^ = R{w, q)/[l — V{q) Qqo], where the 
function R{u), q) is defined as 

+ V{q){Q['^Q\^ -QoQQu)qkqi. (44) 

Notice that only the conventional tensor q\q appears in 
Eq. (44). The anomalous tensor Q^^^ docs not appear 
because it is transverse, as will be discussed in Sec. VA. 
The dispersion relation for the collective modes is deter- 
mined by the equation R{uj, q) = 0. 

In the limit of small \q\, we can simplify Eq. (44) by 
keeping only the nonvanishing terms Qoo and Qki in the 
first and second lines of Eq. (44), whereas the terms Q^.g-' 

and Q^^i vanish in this limit. Using Eq. (40) for 

we find the dispersion relation for the collective modes 

at T = 0, 



V{q 



(45) 



In the case \q\d <C 1, it is appropriate to use the 3D 
limit VzT){q) from Eq. (42) for the Coulomb interaction 

in Eq. (45) and, thus, we find the following spectrum of 
collective modes [48, 49, 56]: 



2 ^^Ib^l + ^lll 
= ^—^ 



= UJ 



qnq 



ah ' 



(46) 



However, this limit does not correspond to Sr2Ru04. 

At higher temperatures T Tc, there may be 
long-wavelength oscillations with the acoustic spectrum. 
These oscillations are neutral and consist of supercur- 
rent oscillations compensated by oscillations of the nor- 
mal current, as discovered by Carlson and Goldman [60]. 
However, at T <C Tc only the plasma mode survives, be- 
cause the normal density is exponentially suppressed. 



IV. CONVENTIONAL NONCHIRAL 
ELECTROMAGNETIC RESPONSE 

A. Transverse and longitudinal tensors 

First, we describe general properties of the tensors Qki 

and Q^^Q- It is convenient to separate the symmetric 

tensor Qki into the longitudinal q\i and transverse Q^; 
parts defined by the following relations 



Qki = qIi + Qii, 

QkQkl — IkQkh Qklll = QklQl, 

QkQki = QkiQi = 0- 



(49) 
(50) 
(51) 



Equations (50)-(51) have the following general solution 



kl 



^ krQj'Qs^ si 
QuQuvQv 



Qki — Q 



kl- 



(52) 



In the dynamic limit, w ^ and q 0, which is relevant 

for optical measurements, the paramagnetic term Q^^K 
defined in Eq. (A2), vanishes, and Qki is given by the 
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diamagnctic term (40). Then, the transverse and longi- 
tudinal parts of the tensor Qfez are 



written in the following form: 



{n-q){q-n) 
qnq 



kl 



go 



lb ( {n-q){q-n) 



qnq 



(53) 

(54) 



5p 

3 
P 



-iq ■ P, 

■ A- iujP, 



(59) 
(60) 
(61) 



kl 



The conventional current-density correlation function 
(5^o\ given by Eq. (A3), is an odd function of q and u). 



Here 6p and j are the induced charge and current densi- 
ties, and P is the polarization vector. In Eq. (60) we use 

the shorthand notation [Q"""]/;; = Q^i [see Eq. (53)]. 
(g-. Equation (61) expresses P in terms of the external 

The term Q)^^> satisfies the following identity for small q: ^^^^^^-^^ ggj^ E''^' = -iqA^ + icjA (Wc remind the reader 

that the electromagnetic potential A^^ in our calculations 
represents the external field.) E'^^^ is connected to the 
total electric field which includes the field created 
by other electrons in the system, by the standard relation 



r)(«) - nil n, 

^Ok — ^klll T|i 



QrQrsQs 



(55) 



To prove Eq. (55), we observe that the tensor structure 

of Qki = JvkVi... and Q^q* = Jvk-.- in Eqs. (39), 
(A2), and (A3) is determined by the electron velocities 

V = dS^p/dp. Since the term Q^^^ vanishes at q = 0, 
the leading-order expansion of the integrand is propor- 
tional to VkViqi, i.e., qJ^'q = / VkViqi . . .. Thus, in the 

long- wavelength limit, wc have oc QkiQi, which leads 
to Eq. (55). Equation (A3) gives the following explicit 

expression for QL for small q and finite w at T = 0: 



^kO 



^lb^(n-q)k 



47r 
1 



^0 



dx 



x^sfx^ 



(56) 



■ 2a; + (w i(5)/Ao 2a; - (w + i^)/Ao ' 



We see that Qj.q-' is proportional to the tensor n multi- 
plied by q, which is consistent with Eq. (55). 



B. Conventional nonchiral electromagnetic 
response 

We now show that if we omit the anomalous chiral term 
(29) in Eqs. (33) and (36), we recover the conventional 
response of a superconductor to the electromagnetic field. 
Using Eqs. (33) and (49) and the identity (55), we obtain 
the following expression for the space-time components 
of the conventional response kernel (36): 



(c) 
kl 



Qii 



kl 
^10 



■ to^Uki, 



K, 



(c) 
00 



-QkHkiqi, 



(57) 



K, 



01 



In the long- wavelength limit assumed here, the transverse 
tensor Qjy is given by Eq. (53), and the longitudinal ten- 
sor kki is defined as follows: 



Kkl 



QknQai QkiQ' 



R 



(58) 



with the function R(LU,q) given by Eq. (44). By using 
Eq. (57), the charge and current responses (43) can be 



E^ 



(62) 



where *s is the tensor of dielectric permeability. If we 
expressed the polarization vector P in Eq. (61) as a lin- 
ear function of the total electric field then the cor- 
responding proportionality tensor would be the dielec- 
tric susceptibility tensor Kki [61]. However, because Eq. 
(61) expresses P as a function of the external field E'^^^, 
the corresponding tensor kki is the RPA-renormalized di- 
electric susceptibility tensor, which includes the effect of 
charge screening via the RPA diagrams shown in Fig. 2. 
The tensor kki in Eq. (58) can be simplified in the limit 

of small q. In this case, the terms involving Q^.^ in the 
numerator and denominator of Eq. (58) can be neglected 

(s) 

relative to the other terms because Qj,q vanishes at qr ^ 
0. Using this approximation and Eq. (44), we rewrite 

Eq. (58) in terms of the longitudinal tensor qJ-^ (54) as 
follows: 



ah 



a{w,q) 



00 



qQq + w^Qoo - V{q) QaoQ-Q-q 



(63) 

(64) 



Here, the dimensionless function a{oj, q) describes screen- 
ing of charge in the static case and plasma oscillations in 
the dynamic case. 

Indeed, in the static limit w = 0, the function Qoo, 
given by Eq. (Al), is proportional to the density of states 

goo(w = 0) = -2e^No/d, No = me/27r, (65) 

where Nq is the 2D density of states per spin. Then 
Eq. (64) for a; = reads 



,(-\co = 0,q) = -—^ 



(66) 



{q-n-q) l + q^j./q'^' 
where qTF is the inverse Thomas-Fermi screening length 

4f = STTe'^No/d. (67) 
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Thus, Eq. (66) describes the electrostatic screening of 
charge in an anisotropic conductor. 

In the dynamic case lo 3> wfIqI, we neglect the first 
term in the denominator of Eq. (64) and obtain 



Sd) 



(68) 



Equation (68) exhibits a resonance, signifying a diver- 
gence of the charge response, when the frequency to ap- 
proaches to the plasma frequency defined in Eq. (46). 

Equations (58)-(61) coincide with the results of 
Ref. [46] [Cf. with Eqs. (28) and (29) of Ref. [46], which 
were obtained using an exact solution for vertex func- 
tions]. Note that the transverse part of the current re- 
sponse in Eq. (60) is not affected by the collective modes 
[45] , whereas the charge response Eq. (59) is strongly af- 
fected by the collective dynamics of the superconducting 
phase and the Coulomb interaction. Equations (59) and 
(60) satisfy the continuity equation lo Sp = q ■ j and are 
invariant with respect to gauge transformations of the 
electromagnetic field. In anisotropic superconductors, it 
is not practical to separate A and E'^^^ into longitudinal 
and transverse components, because such separation does 
not diagonalize the response equations (59) and (60), un- 
like in the isotropic case. 



V. ANOMALOUS CHIRAL 
ELECTROMAGNETIC RESPONSE 

A. Anomalous current-density correlation function 

The chiral anomalous current-density correlation func- 
tion q'^i^q is given by Eq. (A4). By separating the factors 



Qx and Qy, it can be written as 



Q 



(a) 
kO 




(69) 



where r]xy — —rjyx = 1 is the 2D antisymmetric tensor, 
and Qk is the fcth component of the vector 



q — z X q. 



The anomalous term q'^^J (69) is transverse: qkQ 



(70) 



(a) 
fcO 



0, 



unlike the conventional term Q^p* in Eqs. (55) and (56). 

In the limit q 0, we obtain the following expression 
for the function Q{u!) from Eq. (A4): 



e(a;) 



me J (27r)3 4E^ 



l-2f{E)] 



1 



1 



2E + u} + i5 2E~uj~iS 



(71) 



which was found before in Ref. [31]. Changing the vari- 
able of integration to x = E/Aq and taking into account 




FIG. 3: (color online). Frequency dependence of the func- 
tion I{u>) given by Eq. (73). This function determines the 
frequency dependence of the ac Hall conductivity (104) of a 
chiral px + ipy superconductor in the dynamic limit. The 
solid (red) and dashed (blue) lines show the real and imagi- 
nary parts of I{w). 



that the integral converges near the Fermi surface, we 
express Eq. (71) in the following form at T = 0: 



dx 



^/x^~l a;2 - [{iu + i5)/2Ao\ 



(72) 
(73) 



We restored the dimensional constants in Eq. (72). The 
sign of Q is determined by the chirality Sxy of the su- 
perconducting condensate. The function /(w) has the 
following asymptotic behavior [62]: 



1, w<Ao. 

(74) 

The crossover between the two limiting cases occurs at 
LD ~ 2Ao, where the photon energy is equal to the bind- 
ing energy between two electrons in a Cooper pair. The 
real and imaginary parts of /(w) are shown in Fig. 3. 
The imaginary part is zero for w < 2Aq and diverges as 
1 / \/cj~^^"2Ao when w approaches 2Ao from above. The 
real part is approximately constant at low frequencies 
and diverges as l/\/2Ao — oj when u approaches 2Ao 
from below. Both real and imaginary parts go to zero 
at cj — > oo. 



B. Anomalous chiral electromagnetic response 



Now we collect the terms that contain the chiral 
anomalous correlator Q^j^J in the kernel Kfj_i, (36). First, 

we notice that Q^.g'* does not appear in the denominator 
Qapq'^q^ of Eq. (36), because the denominator is lon- 
gitudinal, whereas Q^.g'* is transverse. Then, Eq. (36) 

(a) 



contains only the linear and quadratic terms in 



fcO 
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The quadratic in Q^",^ term may be called the "double- 
anomalous." Although this term has an anomalous ori- 
gin, it is not chiral, i.e., it does not change sign under 
the time-reversal operation or when the chirality Sxy of 
the order parameter changes sign. This term turns out 
to be small in the relativistic parameter {vp/c)'^ and is 
not particularly important. The discussion of this term 
is deferred to Appendix C. 

In the rest of the paper, we concentrate on the linear in 
Qfco^ terms in the response. Using Eq. (36), one obtains 
the following expression for the space-time components 



of the anomalous chiral kernel K^iJ 



(«). 



(a) 
kl 



QkWi - WkQl 
R 



(a) 
00 



(75) 



R 



Here, R, ©, and qk are given by Eqs. (44), (69), and (70), 
and the vector W is defined as 



Wk = Qlqi + QIq-'w. 



(76) 



By using Eq. (75), the anomalous charge and current 
responses (43) can be written in the following form: 

= -iq.p(«)^ (77) 

= iq X M^"> - iLoP'^''\ (78) 

where Af*^") and P^"-' can be identified as the chiral mag- 
netization and polarization [21, 61] and are given by 



p(o) 



-iew 



[q X E"" 



izQ 



u)R 
{W ■ E"''*) 

R ■ 



W 

R ^ 



(79) 
(80) 



Here, B^^^ is the z component of the external magnetic 
field B'^'^*, which is related to E'^^'^ via Maxwell's equa- 
tion q X E"""^ = LoB^'^K Both M(") and P^"'^ in Eqs. (79) 
and (80) are proportional to the chiral response function 
Q (72) and, thus, change sign when the chirality Sxy of 
the order parameter changes sign. The peculiar feature 
of chiral superconductors is that the electric polarization 
(79) is induced by the magnetic field B^, whereas 
the magnetization Af(") (80) couples to the longitudi- 
nal electric field. The first and second terms in Eq. (78) 
give the transverse and longitudinal components of the 
anomalous current, respectively. 

In the limit of small |q|, by using Eqs. (54) and (56) 
and the condition /(O) = 1, the vector W in Eq. (76) can 
be written as 



Wk^qiQi 



I{co)^{n-q)k, 



(81) 



where /(w) is given by Eq. (73). In the same limit, only 
the terms proportional to Qoo survive in the function R 
given by Eq. (44). Using Eq. (Al) at q ^ 0, we find that 



(yoo(w) = J — -((w) 



(82) 



where I{uj) again is given by Eq. (73). Then, with the 
help of Eqs. (81), (82), and (44), we find that 



W a{u!,q)d 



R 



4e2 



me 



[n ■ q), 



(83) 



where the function a(uj,q) is given by Eq. (64). Using 
Eq. (83), the expressions (79) and (80) for P(«) and M(«) 
can be written as 



a{ijj, q)d [qx E'^^'^] ■ z 



(n • q) (84) 



= -izG^h^(g-n •£;-*). (85) 



4e2 



me 



The static and dynamic limits of Eqs. (84) and (85) can 
be obtained by taking appropriate limits (66) and (68) 
for the function a{uj,q). 



C. Discussion of the anomalous magnetization 
current 

By comparing Eq. (85) for the anomalous magneti- 
zation with the expression for the induced conventional 
charge density from Eqs. (59), (61), (63), and (54), one 
can notice that Af^") can be expressed in terms of 5p 



Af(°) 



4me 



(86) 



Then, the first term in Eq. (78), the magnetization cur- 
rent, can be written in real space as 



j^^(w,r) = ^ [z X V(5p(w,r)]. 



(87) 



At a; = 0, Eq. (87) reduces to expression (3) for the 
Mermin-Muzikar current [32], which was discussed in 
Sec. I. However, at high frequencies, the anomalous 
magnetization current (87) is suppressed by the function 
I{uj) oc (Ao/w)2 [see Eq. (74)]. The general reason for 
this suppression is that physical manifestations of the 
low-energy Cooper pairing with the angular momentum 
Lz = h should fade away at high frequencies a; l^- Aq. 
Indeed, the electron states at high energies are essentially 
the same as in a normal metal, so the properties of the 
system should approach those of a normal metal at high 
frequencies. The suppression of the anomalous magneti- 
zation current at high frequencies is one of the reasons 
why the Kerr angle is so small [33]. Equation (87) gives 
an important generalization of the Mermin-Muzikar cur- 
rent to arbitrary frequencies. 

In order to calculate the magnetization current cor- 
rectly, it is very important to take into account the term 
Q^Q^ in Eq. (81). This term vanishes at u; = and can 
be neglected at low frequencies, as was done in Rcfs. [25 - 
29]. However, at high frequencies, Eq. (56) shows that 
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Qak ^ ^Ik/^i and, in the loading order of approxima- 
tion, the second term of the sum in Eq. (81) cancels the 
first term. The remaining difference gives the small fac- 
tor I{uj) cx (Ao/w)^. The term Qoo is also proportional 
to I(uj) in Eq. (82), so the factors in W and R can- 
cel out in Eq. (83), but 6(a;) produces the factors 

in Eqs. (85), (86), and (87). If Qo^ ^erc neglected in Eq. 
(81), then W would have a constant value at high fre- 
quencies, and one would incorrectly conclude that Eq. (3) 
is valid for arbitrarily high frequencies. The asymptotic 
cancellation of the two terms in the sum in Eq. (81) can 
also be shown using the Ward identities [46] . 

D. Discussion of the anomalous polarization 
current 

The second term j'^p' in Eq. (78) is determined by the 

chiral polarization P^°'\ Although j^p^ is formally pro- 
portional to w, it contributes equally to the total current 
even at low frequencies. Indeed, Eqs. (79) and (84) show 
that the chiral polarization diverges at low frequencies as 
l/u>. This divergence exactly cancels w in Eq. (78), which 
results in a finite contribution from the chiral polariza- 
tion to the current even at low frequencies. Both magne- 
tization and polarization currents equally contribute to 
the Hall conductivity tensor, which is discussed below in 
Eqs. (100) and (103). The anomalous polarization cur- 
rent jp ' was often omitted in the previous papers in the 
low frequency limit. 

Equations (77), (79), and (84) show that the magnetic 
field Bz induces an anomalous electric charge, as men- 
tioned in Sec. I. The polarization current j^'p' in Eq. 
(78) is necessary to satisfy the continuity equation for 
the anomalous charge in Eq. (77). 

E. Anomalous chiral efFective action 

Substituting the components (75) of the anomalous 
tensor k'^J into Eq. (35), we obtain the chiral part of 
the gauge-invariant effective action for the electromag- 
netic field as 

Given that Eq. (75) is written for the real frequency w, we 
write Eq. (88) for the real frequency as well. The anoma- 
lous charge and current responses, given by Eqs. (77) 
(80), can be obtained by taking the appropriate varia- 
tions of Eq. (88). The causality of the response function 
should be properly addressed, as discussed in Ref. [50]. 

Unlike the conventional action for the electromagnetic 
field, the anomalous action (88) involves a product of 



the electric and magnetic fields. Thus, the anomalous 
action (105) breaks the TRS, because E ^ E and B 
—B upon the time-reversal operation. Equation (105) is 
manifestly gauge-invariant and is a replacement for the 
Chern-Simons-type term (1) after integration out of the 
superconducting phase <I>. The calculation of this action 
is one of the central results of our paper. A simplified 
alternative derivation of the effective action is also given 
in Appendix C. 

Anomalous effective actions of the forms similar to 
Eq. (88) were obtained for chiral superfluids in Refs. [25, 
26, 28] . Coupling between the electric and magnetic fields 
was discussed in Ref. [63]. However, because these calcu- 
lations were performed in the low-frequency limit, they 
did not obtain the factor /(w), which suppresses the chi- 
ral effects at high frequencies, as discussed in Sec. VC. 
The factor R in the denominator of Eq. (88) represents 
the collective modes of the system. References [25, 26, 28] 
did not take into account the Coulomb interaction, which 
is appropriate for electrically neutral superfluids, such 
as ^He. Thus, the collective modes in the denominator 
of the effective actions in these papers were the acous- 
tic modes of the superconducting phase, which can be 
obtained by setting V{q) = in Eq. (44). However, 
in the case of a charged superconductor, the collective 
modes are the gapped plasmons represented by the func- 
tion R in Eq. (88) (see also Ref. [28]). One should keep 
in mind that, because we have integrated out the inter- 
nal Coulomb potential, the electromagnetic field in Eq. 
(88) is the external one. Depending on the physical con- 
text, it may also be useful to study the chiral response 
with respect to the total electromagnetic field, which in- 
cludes both external and internal fields. This problem is 
discussed in Sec. VI. 



VI. CONDUCTIVITY TENSOR OF A CHIRAL 
SUPERCONDUCTOR 

A. Response to the external versus total electric 
field 

In order to calculate the polar Kerr angle (see Ap- 
pendix D) and other observable experimental quantities, 
one needs to know the conductivity tensor a of a chiral 
superconductor, which is defined by the standard relation 

j = 'ff. E^°\ (89) 

Here, E*'°^ is the total electric field inside the supercon- 
ductor. In principle, the tensor can be extracted from 
Eqs. (60) and (78) for the conventional and chiral cur- 
rents. However, the polarization and magnetization vec- 
tors in these equations arc expressed by Eqs. (61), (79), 
and (80) in terms of the external electric field E'^^'^ . These 
equations give the linear response relation in the form 

j = d- E^"""- (90) 



13 



with a different tensor a [42], which inchides the renor- 
malization due to the RPA diagrams shown in Fig. 2. 

The linear response relation in the form (90) is phys- 
icahy transparent, because it gives a direct response of 
the system to the external perturbation. However, in 
general, the response of the system to the external field 
depends on the geometry of the sample and the exper- 
imental apparatus. Only in the idealized case of an in- 
finite uniform system can the problem be solved by the 
Fourier transform. To deal with this problem, the stan- 
dard approach in the electrodynamics of continuous me- 
dia [61] is to use the constituency relation (89), which 
expresses the response of the media to the total electro- 
magnetic field in terms of the conductivity tensor char- 
acterizing the material. The constituency relations are 
substituted into Maxwell's equations as the source terms, 
and then Maxwell's equations are solved with the bound- 
ary conditions appropriate for the experimental setup. 
This is how, for example, one can take into account the 
Meissncr screening, which was not included in the RPA- 
renormalized linear response given by Eq. (90). 

In order to obtain the proper constituency relations, 
we need to transform the results of the paper to the 
form (89). This can easily be done by noticing that, 
in Eq. (12), the internal (induced) electric potential 
appears next to the external potential A^^. Thus, the 
combined electric potential corresponds to the total elec- 
tric field E^°^ inside the superconductor. Then, Eq. (23) 
gives the effective action for the total electromagnetic 
field. By taking a variational derivative with respect to 
the total electromagnetic field, one arrives at Eq. (89). 
Formally, this means that we do not integrate out the 
internal field ip and do not perform the RPA renormal- 
ization of the response kernels Qafj in Eq. (33). Thus, 
the transition between Eqs. (90) and (89) can be accom- 
plished by setting V{q) = and replacing E""^^ with E^°^ 
in the electromagnetic response. The Coulomb poten- 
tial V{q) appears explicitly only in the function R{uj, q) 
defined in Eq. (44). Therefore, one should replace this 
function with the bare one [64] , 

R{uj,q) = Qkmqi + '2u;Q^ok Ik + '^'^Q 00- (91) 

Using this prescription and Eqs. (60) and (78), one can 
easily obtain the conductivity tensor 



(92) 



which consists of the conventional (nonchiral) ct^"^^ and 
the anomalous (chiral) a^"^ contributions, as discussed 
below. 



of the conductivity tensor as follows: 



?(c) 



±(2) 
'^kl 



ICO 



-lUl- 



(QkoQoi ~ QkiQoo) 



R 



(93) 
(94) 

(95) 



0, the terms Q^^fj 



where the function R{uj,q) is given by Eq. (91). 

In the dynamic limit w 7^ and q 
in the numerator of Eq. (95) and the terms proportional 
to q in Eq. (91) vanish. As a result, the expression for 
the conductivity tensor is simplified 



'kl 



?(2) _ 



1 



'kl 



kl- 



(96) 



From Eq. (96), we find the total conductivity tensor in 
the long wavelength limit as follows: 



(c)__QiP 



^ki = - 



<b 
ojI, 




(97) 



1(1) 



where we used Eq. (40) for Q^; . 

Combining Eq. (97) with the standard formula for the 
dielectric permeability tensor *e , 



e = 1 + CT, 



(98) 



we obtain 



e(u;) 




(99) 



Equations (97) and (99) represent the standard Drude 
response of an anisotropic conductor. 



C. Anomalous chiral conductivity tensor 

The anomalous chiral response of the system to the 
total electric field is obtained from Eqs. (77)-(80) by re- 
placing R ^ R, where the function R{uj,q) is given by 
Eq. (91). From Eqs. (78)-(80) for the current response, 
we obtain the anomalous chiral conductivity tensor 



QkWi - WkQi 



(100) 



B. Conventional nonchiral conductivity tensor 

Following the prescription of Sec. VIA and using 
Eqs. (58), (60) and (61), we obtain the conventional part 



where q and Wk are given by Eqs. (70) and (76). Notice 
that the anomalous conductivity tensor is antisymmet- 



'kl 



'^ik'il) ^^"^ represents the intrinsic Hall 



conductivity of a chiral superconductor. 
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Using Eq. (81), we can rewrite the anomalous conduc- 
tivity tensor (100) as 

Cqxqz I ,(102) 





Qqyqz -QqxQ. 



where C = cj^/cj^^ is a small parameter representing 
anisotropy of the band dispersion in a Q2D metal. The 
tensor components and in Eq. (102) are small, 
unless ql > Q'^{ql+ql). 

Now we concentrate on the a'^xy component of the 
anomalous conductivity tensor (102) 



ctW = -6(0;) 



4TrR{u), q) 



q\=ql+ql. (103) 



It was argued in Refs. [31, 33] that the Hall conductiv- 
ity of a chiral superconductor can be obtained by omit- 
ting the last term in the brackets of Eq. (4), which gives 

<y^y = 6(w). However, Eq. (103) shows that the system- 
atic integration out of the superconducting phase $ pro- 
duces a different expression for a^xy , which differs from 
0(a;) by the additional factor proportional to This 
factor involves the in-plane wave vector of the electro- 
magnetic field and, in the appropriate limit, makes a^y 
vanish when ^ 0. This fact reflects the cancellation 
of the Hall effect [34] for a uniform (in the plane) sys- 
tem discussed in Sec. I. Equation (103) is the central 
result of our paper and will be used in Sec. VII to es- 
timate the observable Kerr effect. At low frequencies 
ti) <C Ao, where / w 1, Eq. (103) agrees with the cor- 
responding results of Refs. [25, 28]. However, at higher 
frequencies a; > Aq, the function I{uj) exhibits the non- 
trivial frequency dependence shown in Fig. 3. This be- 
havior originates from the tendency to cancel between 
the current- current Qki and current-density Q^^q polar- 
ization functions in Eq. (81). Thus, it is essential to take 

Qfeo* iiito account at high frequencies (see discussion in 
Sec. VC). 

Let us discuss Eq. (103) in different limits. First, we 
consider the static limit w = 0, set = 0, and then 
take q\\ 0. Using Eqs. (72) and (91), we see that 

cancels out, and we find that ai^ = e'^/2hd, which 
is reminiscent of the quantum Hall effect [14, 31, 33]. 
However, this formally calculated value does not corre- 
spond to an observable dc Hall effect. The static limit de- 
scribes the system in thermodynamic equilibrium, where 
an applied electric field causes an inhomogeneous equi- 
librium redistribution of the electron density. As a result 
of the nonzero Cooper-pair angular momentum, the in- 
homogeneous electron density produces the equilibrium 
magnetization current Eq. (87). However, the total mag- 
netization current flowing through any cross section of 



the sample, including the bulk and the edges, is zero, 
because the current is solenoidal. Therefore, the total 
current measured by an ammeter is zero. Thus, the for- 
mally calculated does not represent a measurable 
Hall effect [65]. 

The experimentally relevant limit is the dynamic limit 
with w 7^ and small q. Taking this limit in Eq. (103) 
and using Eqs. (91) and (82), we find 



xy 



2^2 ' 



(104) 



where vp is the in-plane Fermi velocity. Equation (104) 
gives the ac Hall conductivity, where the frequency de- 
pendence of 8(ti>) is given by Eqs. (72) and (73). Equa- 
tion (104) differs from the expression ai'ij = Q{oj) ob- 
tained in the previous papers [31, 33] by the small factor 



D. Anomalous effective action and charge response 

Equation (88) gives the anomalous effective action as a 
function of the external electromagnetic field. Transfor- 
mation of this action to the total field is straightforward 
by replacing R [64], 



qia) 



R{q) 

q-n-E^°\q)]Bl°\-q) 
- vl.{q • n • q)/2 



(105) 



The difference between Eqs. (88) and (105) is that the 
former involves the gapped plasmon modes represented 
by R (44) in the denominator, whereas the latter involves 
the acoustic modes of the siiperconducting phase repre- 
sented by R (91). This difference occurs because the two 
effective actions are written using the screened and un- 
screened electric fields, whereas the magnetic field is not 
screened by the internal Coulomb potential. For small q 
and high u), this difference amounts to using R Ri Qoquj^ 
vs R ^ Qoo[lu'^ — u;2j^(<7 ■n-q)/q'^], where the last term 
represents the momentum-dependent plasma frequency. 
The effective action is also discussed in Appendix C. 

By taking a variational derivative of the action (105) 
with respect to the total field ^o°*i '^^ ^^'^ the anomalous 
electric charge induced by the magnetic field [66], 

Spi-) = -!^^^q.n.q)Bl°\ (106) 
In the static limit, Eq. (106) reads 

c2 



dp 



(a) 



2hdc 



(107) 



Equation (107) was derived earlier in Refs. [27, 29] and 
was mentioned in Sec. I. It is reminiscent of the Streda 
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formula for the quantum Hall effect [67]. However, at 
high frequencies, Eq. (106) shows that the induced elec- 
tric charge is significantly reduced relative to the static 
limit (107) 



2hdc 



v%{q-n-q) 

2u;2 



Bl 



(108) 



It is instructive to calculate how much electric charge 
would be induced by the static magnetic field of one su- 
perconducting vortex. If we take an integral / dx dy of 
Eq. (107), the left-hand side gives us the total induced 
electric charge, and the right-hand side gives the total 
magnetic flux. Taking the latter to be one supercon- 
ducting flux quantum 0o = hc/2e, we find that the total 
induced electric charge in a vortex JQvortex is 



*^^vortex — ^xy 



(109) 



per layer, i.e., the superconducting vortex has a frac- 
tional electric charge. Equation (109) was obtained by 
Goryo [27], who also pointed out that a vortex has a frac- 
tional angular momentum. However, the anomalous in- 
duced charge (109) would be screened by the conventional 
screening mechanism [66] , and its experimental measure- 
ment may be challenging. It should be emphasized that 
the charge density (107) is not concentrated in the vortex 
core at the coherence length, but extends to the London 
penetration length, where the magnetic field is present. 
A possible experiment for detection of the anomalous in- 
duced electric charge is discussed in Sec. VII B. 



VII. EXPERIMENTAL IMPLICATIONS 



A. Polar Kerr effect experiment 

In this section, we apply the derived theoretical results 
to the interpretation of the Kerr effect measurements [3] . 
We use the relationship between the Hall conductivity 
(Jxy and the Kerr angle 9k as presented in Appendix D. 

First, we briefly discuss the estimates of the polar Kerr 
angle in the previous literature [31, 33], where the for- 
mula a^xy = 0(w) was used without the additional fac- 
tor appearing in Eq. (104). To estimate the Kerr angle, 
Ref. [31] implicitly assumed that the refraction coefficient 
n is real by using the value n(n^ — 1) = 3 quoted in Ref. [3] 
and expressed Ok in terms of aij.y via Eq. (D2) . The the- 
oretical estimate given in Ref. [31] is 9k ~ 230 nrad, 
which is of the same order of magnitude as the experi- 
mental value of 65 nrad [3] . In Ref. [33] , a more detailed 
estimate was presented, using Eq. (D4) for the refrac- 
tion coefficient n and discussing different limits u) > Wp 
and u> < Up. However, the factor Coo was overlooked in 
Ref. [33] as well as the presence of the real part <T'^y{uj), 
which gives the primary contribution to 9k for w < u)p { 
see Appendix D). Overall, the numerical estimate of 6k 
given in Ref. [33] is of the same order of magnitude as in 
Ref. [31]. 



However, as discussed in Sec. I, Refs. [31] and [33] 
did not take into account the self-consistent dynamics 
of the superconducting phase $ and, thus, missed the 
additional factor in Eq. (104), which makes the Hall con- 
ductivity dependent on the wave vector . In this case, 
strictly speaking, one cannot use the equations for 9k 
presented in Appendix D, because they were derived for 
the normally-incident infinite plane wave with = 0. 
However, as shown in Fig. 1, the experiment [3] was per- 
formed with a tightly focused Gaussian laser beam, which 
has nonzero Fourier components with g|| ^ 0. Solving 
the boundary value problem for a reflection of a flnite- 
size Gaussian beam from a chiral superconductor and 
determining the Kerr angle for polarization rotation is 
a complicated problem, that is beyond the scope of this 
paper. Nevertheless, to make a crude estimate of the 
Kerr angle, we can use the equations from Appendix D 
and the Hall conductivity given by Eq. (104) with the 
replacement — > l/l, where I is the typical transverse 
size of the Gaussian beam ( see Fig. 1). Taking into ac- 
count that LV = c\q\, where the wave vector \q\ = 27r/A 
is related to the wavelength of light A, we can roughly 
estimate the Hall conductivity in Eq. (104) as 



V c 



(110) 



Using the values vp = 5.5 x 10^ m/s for the 7 sheet of the 
Fermi surface [10], c = 3 x 10^ m/s, A = 1.55 ^m, I ~ 25 
Hm [3], and Aq = 0.8 mcV [20] in Eqs. (D7) and (110), 
we estimate the Kerr angle to be 0k ~ 10"^'' rad. This 
estimate for 6k is about 6 orders of magnitude smaller 
than the experimental value of 65 nrad. The strong sup- 
pression of the Kerr angle relative to the previous esti- 
mate originates primarily from the small relativistic fac- 
tor (wf/c)^ in Eq. (110). 

In the rest of this section, we speculate about possi- 
ble ways of resolving the discrepancy between the theory 
and experiment. As discussed in Sec. I, the general ar- 
guments of Ref. [43] show that the intrinsic Hall conduc- 
tivity should vanish for a spatially homogeneous uniform 
system. One way to break the translational symmetry is 
by taking into account the finite size I of the laser beam. 
Then, inevitably, the Kerr angle acquires a dependence 
on I as shown in Eq. (5), which is a very robust the- 
oretical result. This proportionality relation should be 
checked experimentally. 

Although the dependence on I introduces a small fac- 
tor, there may be mechanisms for enhancement of the 
response of the system, which may compensate for this 
additional smallness. One possibility is a resonance with 
the plasma modes. Equation (87) shows that the mag- 
netization current is produced by the gradients of elec- 
tron density. Therefore, the problem of the chiral cur- 
rent calculation reduces to how much electron charge is 
induced on the surface of a crystal by the inhomogeneous 
laser beam. The continuity of the electric field lines in a 
Gaussian beam in vacuum requires that the electric field 
components must be present at the periphery of the 



beam, even though the electric field is nominally polar- 
ized along at the center of the beam. This is well 
known for the fiber-optics modes, which have similari- 
ties with the Gaussian mode in vacuum [68]. When the 
Gaussian beam hits the sample (see Fig. 1), the Ez com- 
ponents induce electric charges of opposite signs on the 
two sides of the beam in the x direction, which induce 
the electric current jy oc QxE^ according to Eq. (87). The 
induced electric current jy generates a magnetic field 
and a reflected electromagnetic wave with the Ey polar- 
ization. Taking into account that E^ oc {q\\/\q\)Eg;, we 
obtain the same result as in Eq. (104). However, because 
the external electric field E^ directly couples to the in- 
duced electron charge, the response of the system would 
be enhanced by the factor a^'''(u;,q) in Eq. (68) when 
the frequency of light lo is in resonance with one of the 
plasma modes. As discussed in Appendix D, the exper- 
imental frequency to is in between the upper and lower 
plasma frequencies uJab > to > Wc and, thus, to can be in 
resonance with one of the plasma modes (46) for some 
vector q. Besides, because of the boundary at z = be- 
tween the sample and vacuum, it may be necessary to 
solve for the plasma modes more accurately by taking 
into account the surface plasmons as well. 



Another reason for the discrepancy in the magnitude 
of the Kerr effect may be the idealization of our theo- 
retical model. Indeed, we considered the electromagnetic 
response of a clean superconductor in the presence of 
particle-hole symmetry. Within this model, Eqs. (Al)- 
(A4) for the response kernels are written assuming mo- 
mentum conservation for the electrons. A more realistic 
model has to take into account the effect of impurities, 
in which case the momentum conservation does not hold. 
Indeed, one of the reasons for smallness of 9k within 
our model can be traced back to the cancellation of the 
two terms in Eq. (81) at high frequencies. As a result, 
the right-hand side of Eq. (81) is small in the parame- 
ter (Ao/huj)'^ ^ 10~^. If we take into account impurities, 
the complete cancellation in Eq. (81) might not hold, and 
the estimate for 6k would be significantly larger. Indeed, 
in the derivation of this equation, we assumed that the 
paramagnetic kernel Q^/ given in Eq. (A2) vanishes in 
the dynamic limit [lo and q — > 0). However, it is 
well known from the classic paper by Mattis and Bardeen 
[69] that Qj^i does not vanish when impurity scattering 
is taken into account. This term plays a significant role 
by taking away some spectral weight from the London- 
Meissner supercurrent, which makes the superconduct- 
ing gap visible in optical measurements. For Sr2Ru04, 
this fact was experimentally confirmed in Ref . [70] , which 
found that the superfluid density at T = is reduced 
by 22%. Therefore, the cancellation in Eq. (81) in the 
presence of disorder should be investigated in a future 
theoretical work. 
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FIG. 4: (color online). Schematic setup of a proposed exper- 
iment to measure the electric charge induced by a magnetic 
field perpendicular to the layers in Sr2Ru04, as given by Eq. 
(107). A current in the solenoid or the coil produces a mag- 
netic field Bz inducing an electric charge at the center, which 
is carried by the radial currents shown by the arrows. The 
induced charge is detected by the AFM tip. 



B. Proposed experiments 

As discussed in Sec. VI D, the magnetic field compo- 
nent Bz applied perpendicular to the conducting planes 
of Sr2Ru04 is expected to induce an electric charge. The 
effect is stronger in the static limit than in the dynamic 
limit. In this section, we propose a conceptual exper- 
iment to verify this effect. The schematic experimen- 
tal setup is shown in Fig. 4. A miniature solenoid or a 
coil carrying electric current creates a magnetic field with 
the Bz component. The current and the magnetic field 
are modulated at a low frequency (a few kilohertz), so 
that the quasistatic formula (107) is applicable [71], but 
the lock-in measurements are enabled. The slowly alter- 
nating magnetic field Bz induces an alternating electric 
charge near the center of the coil, which is produced by 
the radial alternating currents, shown by the arrows in 
Fig. 4. The induced electric charge is detected by a sen- 
sitive atomic force microscope (AFM) tip, shown as the 
pointed vertical object in Fig. 4. The AFM tip should 
be made from a nonmagnetic but highly electrical polar- 
izable material to reduce the response to magnetic fields 
and enhance the response to the electric charge. To in- 
crease sensitivity, measurements should be performed at 
the lock-in frequency of the current in the coil. The dif- 
ference of the signals above and below the superconduct- 
ing transition temperature should be taken in order to 
subtract the background effects. A noticeable increase 
of the signal below Tc would signify the observation of 
the effect. Then, the measurements can be performed 
at different locations on the sample. If a superconduct- 
ing domain of the opposite chirality is found, the effect 
would change sign, i.e., the sign of the induced electric 
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charge would change vs the sign of the apphcd magnetic 
field. It may be possible to map out the chiral domain 
boundaries in this way, assuming that the domains are 
sufficiently static and do not fluctuate in time too fast. 

VIII. CONCLUSIONS 

In this paper, we studied the electromagnetic re- 
sponse of a Q2D chiral px + ipy superconductor, such 
as Sr2Ru04. By integrating out the superconducting 
phase $, we obtained the gauge-invariant effective ac- 
tion of the system in an electromagnetic field. Besides 
the well-known conventional terms, this action also con- 
tains an anomalous chiral term, which breaks the time- 
reversal symmetry. Instead of the Chcrn-Simons-type 
term (1), discussed in the earlier literature [14, 30, 31], 
the gauge-invariant anomalous effective action (105) con- 
tains a product of the electric and magnetic fields. This 
is rather unusual situation, and the modified Chern- 
Simons-type term (105) leads to a number of interesting 
effects. 

By taking variations of the anomalous effective action 
with respect to the electromagnetic potentials Af^, we 
obtained the anomalous charge and current responses of 
the system. We found that the transverse chiral response 
couples to the collective plasma modes, while the anoma- 
lous charge response manifests itself as the electric charge 
induced by the magnetic field perpendicular to the 
conducting layers of Sr2Ru04. In Sec. VII B, we proposed 
an experiment for the detection of this charge. The effect 
is stronger in the static case than at high frequencies. 

The anomalous current response manifests itself as the 
intrinsic Hall effect, i.e., the existence of the antisymmet- 
ric Hall conductivity axy = —(Jyx in the absence of an 
external magnetic field. In the static case, it is a formal 
consequence of the magnetization current J = V x M 
(the Mermin-Muzikar current [32]) in the presence of an 
inhomogeneous electron density induced by an electric 
field. Because the Cooper pairing takes place with the 
angular momentum = ±?i, the variation of electron 
density causes a variation of magnetization. However, 
the equilibrium magnetization currents, obtained in the 
static limit, are solenoidal and cannot be directly mea- 
sured in the transport experiments. In the dynamic limit, 
our calculations show that the magnetization current is 
greatly reduced by the parameter (Ao/w)^. In addition 
to the magnetization current, there is also an anomalous 
polarization current, which originates from the motion 
of the electric charges induced by the magnetic field B^. 
Both the magnetization and polarization currents are of 
the same order and are suppressed at high frequencies. 



We found that the calculated Hall conductivity is pro- 
portional to where is the in-plane wave vector of 
the electromagnetic wave. This is an interesting result, 
which follows from the cancelation of the electric field by 
the supercurrent. The latter was not taken into account 
in the previous papers [31, 33]. As a result, we found 
that the Kerr angle 6k should depend on the transverse 
size I of the laser beam (see Fig. 1) according to Eq. 
(5). Within the considered model here (a clean super- 
conductor in the presence of particle- hole symmetry), we 
found that the overall magnitude of the Kerr angle is 
much smaller than experimentally observed in Rcf. [3]. 
We pointed out two possible reasons for this discrepancy: 
coupling to the plasma resonances and the effect of impu- 
rity scattering. In the presence of impurities, the quasi- 
particle momentum is not conserved, and the results of 
our calculations would change significantly. We conclude 
that, in order to understand the experiment [3], it is nec- 
essary to take into account the effect of disorder. 

Note added. Recently, a preprint by Roy and Kallin 
(Ref. [72]) appeared, in which similar results were found 
using a somewhat different approach. 
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APPENDIX A: POLARIZATION FUNCTIONS 

In this appendix, we present analytical expressions for 
the correlation functions given by Eqs. (24)-(29) and gen- 
eralized to the 3D case as explained in Sec. II F. After 
calculating the Matsubara sums and performing the ana- 
lytical continuation to the real frequency, we find the fol- 
lowing expressions for the finite-temperature causal po- 
larization functions 
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The notation is explained after Eq. (24), and the integra- 
tion over pz runs from ^n/d to n/d. Hero, f{E) is the 

Fermi distribution function. The diamagnetic term q'^j' 
is given by Eq. (39). 

APPENDIX B: COLLECTIVE MODES 

An alternative way of deriving the collective modes [40] 
is to set = in Eq. (23) and study the dynamics of 
the internal scalar potential (f{q) coupled with the super- 
conducting phase ^{q). After an analytical continuation 
to the real frequencies, the effective action becomes 

Ses{^,'p)= (Bl) 

=ii:(^(.)..(.))(£ £)(-;!-,?')■ 

Here, the matrix elements are 

Mil ^ Qkiqkqi + Qoo'^^ + iQok + Qko)ijJqk, (B2) 
M22 = V{q)-^-QoQ, Mi2 = -M2i = - [QoooJ+Qkoqk] ■ 



Dispersion of this mode is determined by setting the de- 
terminant of the matrix M to zero, 



Det(M) = M11M22 - M21M12 = 0. (B3) 



This equation is the same as the equation .R = with R 

given by Eq. (44). 

APPENDIX C: ALTERNATIVE DERIVATION OF 
THE EFFECTIVE ACTION 



In this appendix, we give a simplified alternative 

derivation of the effective action for a chiral supercon- 
ductor after elimination of the superconducting phase 
This derivation is more transparent and easier to 
compare with the effective actions discussed in Refs. 
[14, 25, 26, 28, 29]. 

The starting point is the effective action (23), which is 
explicitly written below in the tensor components 



Here, in the spirit of Sec. VIA, we take the electromag- netic potentials to represent the total electromagnetic 
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field, without separating into the external and internal 
parts, so we drop the terms with the internal scalar po- 
tential from Eq. (23). To shorten notation, it is implied 
in each term of Eq. (CI) that the first dynamical variable 
has the argument q, e.g., Ao(<7), and the second variable 
has the argument — g, e.g., Ao(— q), as in Eq. (23). The 
effective action (CI) is written using the real frequency 
u), as discussed in Sec. VE. 

The superconducting phase $ in Eq. (CI) is the dy- 
namical variable, which should be eliminated by mini- 
mizing S with respect to $ and integrating it out. Be- 
fore doing so, we shift the variable $ by introducing a 
different variable A, 



2e 

^{q) = (l){q) + ^AQ{q). 



(C2) 



Substituting Eq. (C2) into Eq. (CI) and sorting out the 
obtained terms, we find 



kl o— 



(2e)2 



(t> Eh „ (t> „ 
2e LO 2e 



(C3) 

where Ek and arc the electric and magnetic fields. 
The first term in Eq. (C3) corresponds to the Drude re- 
sponse of a metal at high frequencies, as in Eq. (97). 
The second term in Eq. (C3) represents the effective ac- 
tion for the collective variable 0, where the function R 
is given by Eq. (91). The last two terms in Eq. (C3) 
represent the interaction between the collective variable 
(j) and the electromagnetic field. The third term is the 
conventional one, where the function Wk is given by Eq. 
(76). The last term corresponds to the anomalous Chern- 
Simons-type term (1) with the function O being given by 
Eq. (72). Equation (C3) is manifestly gauge-invariant, 
because it contains only the electric and magnetic fields, 
rather than the scalar and vector potentials. 

The cfFcctive action (C3) is quadratic with respect to 
(j). After the elimination of (j), we obtain the final result 

1^/ WuWi\ Ek{q)Ei{-q) 

_ 2i^E^{q)B,{-q) - ^B^B^i-q). 



R 



R 



The first term in Eq. (C4) represents the conventional 

contribution to the effective action discussed in Sec. VI B. 
The second term is the chiral anomalous term (105) 
representing the modification of the Chern-Simons-type 
term due to the dynamics of the collective variable (j). The 
third term in Eq. (C4) is the double- anomalous nonchi- 
ral term, which was briefly mentioned at the beginning 
of Sec. VB. These terms are discussed in more detail 
below. 

The effective actions similar to Eq. (C4) were obtained 
for chiral superconductors after integration out of the su- 
perconducting phase in Refs. [25, 26, 28]. However, the 
focus in these papers was on the low- frequency limit, and 
some terms in the effective action were omitted in this 



limit. The high-frequency limit of the effective action, 
which is relevant for the experiment [3] , was not discussed 
in literature before, except for Ref. [31]. It is remarkable 
that, in the high-frequency limit, all three different func- 
tions G, W, and R « QoqijJ^ are proportional to the same 
function I{co), 

eicv), Wiij), Qoo(w) (X /(w) cx ^, (C5) 

as follows from Eqs. (72), (81), (82), and (74). As a 
consequence of Eq. (C5), the last terms in Eqs. (C4) are 
proportional to and tend to zero as (Aq/w)^ at high 
frequencies. Only the very flrst term in Eq. (C4) survives, 
which represents the Drude response of free electrons at 
high frequencies. The suppression of the last terms in 
Eq. (C4) at high frequencies should be expected, because 
they represent the effect of the chiral superconducting 
state. In the limit Aq/w ^ 0, the superconducting effects 
should vanish, because the low-energy gap Aq cannot 
affect the high-energy behavior, as discussed in Sec. VC. 

By using Eq. (C5), the last double- anomalous term in 
Eq. (C4) can be written as 



<^{aa) 



27r rfrripC^ 



B,{q)B,{-q) 
^ 8tt 



(C6) 



where we restored the dimensional units. Compar- 
ing with the action of a free magnetic field in vacuum 
B{q)-B{—q)/%-K, we see that Eq. (C6) gives an orbital 
paramagnetic contribution to the magnetic susceptibility 
of the system. This contribution is small because of the 
relativistic factor e^/drn^c'^ <C 1 in Eq. (C6). Given that 
this term is nonchiral and small, we do not pay much 
attention to the double-anomalous term in the paper. 

The first term in Eq. (C4), which involves the product 
EkEi, represents the dielectric susceptibility of the sys- 
tem. For small q, the term proportional to WkWi can be 
omitted. The remaining term proportional to Qm, when 
combined with the action of a free electric field in vac- 
uum — ^qE{q) ■ E{—q)/8'jT, gives the the conventional 
effective action for the electric field in the medium 



S^^ = -^Y.E{q).^.E{-q), 



(C7) 



where the dielectric permeability tensor *e is given by 
Eq. (99). If we substitute the internal electric field in the 
form i^'"' = iq(f into Eq. (C7), then we can obtain the 
equation q ■ *e ■ q = for the spectrum of the collective 
plasma modes, which is the same as Eq. (46). 

The effective action (C4) is manifestly gauge invariant, 
and the charge and current response functions, discussed 
in the rest of the paper, can be obtained by taking a 
variation of this action. The function R (91) in the de- 
nominators of Eq. (C4) has zeros at the frequencies of the 
collective modes of $ [64]. If it is desirable to directly 
include the internal Coulomb interaction in the effective 
action (rather than by using Maxwell's equations), the 
kernels Q should be replaced by the kernels Q (33) in 
Eq. (CI) and R should be replaced by R in Eq. (C4). 
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APPENDIX D: RELATIONSHIP BETWEEN THE 
POLAR KERR ANGLE AND AC HALL 
CONDUCTIVITY 

In this appendix, we obtain the equations expressing 
the polar Kerr angle in terms of the ac Hall conductiv- 
ity (see also Ref. [33]). To simplify the notation, we drop 
the index (a) from criy''. 

Let us consider the normal reflection of a linearly po- 
larized electromagnetic plane wave, propagating in the 
z direction, from the (a:, y) surface of a Q2D chiral su- 
perconductor. The reflection coefficient |r| and the polar 
Kerr angle Ok are given by the following equations [21] 



(Dl) 
(D2) 



\n- 



1 



Ox 
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w n [n^ — 1) 



where n is the complex refraction coefficient. The refrac- 
tion coefficient is obtained from the dielectric suscepti- 
bility tensor, which has the general form 



(D3) 



Here, the conductivity tensor ct is given by Eq. (97), and 
'e*^""^ is the background dielectric tensor, which originates 
from the polarizability of the other, nonconducting bands 
in the material. For a plane wave with the electric fleld 
polarized in the x direction, the appropriate component 
of the conductivity tensor is Oxx- Then, the correspond- 
ing refraction coefflcient in Eqs. (Dl) and (D2) is 



ah 



2 ' 



e(-) > 1, (D4) 



where we introduced a shorthand notation eoo ■ Equation 
(97), derived for an ideal clean superconductor, gives only 
the reactive imaginary part cr"^, of conductivity and does 
not contain the dissipative, real part (j^^.. To keep the 
presentation simple, here we discuss only the ideal case 
with cr^^ = 0, but the consideration can be generalized 
to include a'^^ ^ 0. Notice also that Eq. (D2) was de- 
rived assuming that the off-diagonal component Oxy of 
the conductivity tensor is much smaller than the diago- 
nal component cr^x- 

The refraction coefficient n in Eq. (D4) vanishes at 
u) = ujp, where ojp is the so-called plasma edge frequency 

= u>ab/\/e^- (D5) 

Equations (Dl), (D2), and (D4) have different forms for 
the frequencies above and below the plasma edge. For 



uj > ojp, the refraction coefficient n in Eq. (D4) is real, 
so the reflection coefficient (Dl) is < |r| < 1, i.e., the 
electromagnetic wave is partially reflected and partially 
transmitted through the crystal. The polar Kerr angle 
(D2) is given by the following expression in this case 



UJ > LO' 



(D6) 

For frequencies co < cjp below the plasma edge, the re- 
fraction coefficient n (D4) is imaginary. In this case, the 
electromagnetic wave is completely reflected (jrj = 1), 
and the polar Kerr angle (D2) is 



47ra;2 a' 



■>K = 



(D7) 

Equations (D6) and (D7) show that the Kerr angle is 
determined by the imaginary part a'!^^ for lo > LOp and 
by the real part o'^y for < cjp. This statement be- 
comes approximate when the dissipative component a'^^ 
is taken into account, and both the real and imaginary 
parts of a,j.y start to contribute simultaneously to . In 
our theory, a'^y and a'!^y are proportional to the real and 
imaginary parts of the function /(w) deflned in Eqs. (72) 
and (73). The plots of the real and imaginary parts of 
/(w) are shown in Fig. 3, and their asymptotic expres- 
sions are given by Eq. (74) [62]. 

One can notice that the denominators of Eqs. (D6) and 
(D7) vanish at certain frequencies, providing resonance 
enhancement of the Kerr angle Ok- At the plasma edge 
frequency cj = Cijp, the square roots in the denominators 
of Eqs. (D6) and (D7) vanish. At the frequency lo = 
'^ab/V^oo — 1 > Wp, the denominator of (D6) vanishes, 
and Qk changes sign. At this frequency, the reflection 
coefficient (Dl) vanishes because n = 1 in Eq. (D4). Of 
course, in the presence of a^^., these singularities will be 
smeared out. 

The optical properties of Sr2Ru04 in the normal state 
were measured in Refs. [58, 59]; however, the main inter- 
est was in the electric fleld polarization along the c axis. 
The plasma frequencies ^Oah = 4.5 eV and = 0.32 eV 
were obtained from the flts of the data, as well as the re- 
laxation rate 7 ~ 0.5 eV. The value eoo = 10 was quoted, 
but for E'^c. In the Kerr effect measurements [3], the 
frequency of the incoming light was a; = 0.8 eV. This 
frequency is clearly below the in-plane plasma frequency 
u)ab and probably lower than Wp, even for Cqo = 10. In 
this limit, the Kerr angle is given by Eq. (D7). However, 
the effect of the relaxation rate 7, which is comparable 
with w, should also be taken into account. 
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